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Cn . ABSTRACT 

Aims. We present an in-depth study of the two components of the binary system LZ Cep in order to constrain the effects of binarity 
on the evolution of massive stars. 

, Methods. We use a set of high-resolution, high signal-to-noise ratio optical spectra obtained over the orbital period of the system 

to perform a spectroscopic disentangling and derive an orbital solution. We subsequently determine the stellar properties of each 
■ component through an analysis with the CMFGEN atmosphere code. Finally, with the derived stellar parameters, we model the 

*-^H' Hipparcos photometric light curve using the program NIGHTFALL to obtain the inclination and the real stellar masses. 

Results. LZCep is a 9III+ON9.7V binary. It is as a semi-detached system in which either the primary or the secondary star 
j.^ ■ almost fills up its Roche lobe. The dynamical masses are about 16.0 M (primary) and 6.5 M Q (secondary). The latter is lower 

than the typical mass of late-type O stars. The secondary component is chemically more evolved than the primary (which barely 
shows any sign of CNO processing), with strong helium and nitrogen enhancements as well as carbon and oxygen depletions. These 

1 1 _ properties (surface abundances and mass) are typical of Wolf-Rayet stars, although the spectral type is ON 9.7V. The luminosity of the 

secondary is consistent with that of core He-burning objects. The preferred, tentative evolutionary scenario to explain the observed 
t-H ■ properties involves mass transfer from the secondary - which was initially more massive- towards the primary. The secondary is now 

almost a core He-burning object, probably with only a thin envelope of H-rich and CNO processed material. A very inefficient mass 
transfer is necessary to explain the chemical appearance of the primary. Alternative scenarios are discussed but they suffer from more 
\Q . uncertainties. 

, Key words. Stars: early-type - Stars: binarity - Stars: fundamental parameters - Stars: winds, outflows 

' 1. Introduction properties of massive binary systems with well-constrained or- 

y—i .... n . . ,, . bital parameters. 
. . The binary fraction of massive stars is not well constrained, 

> but studies of y oung open cluste rs indicate that it can easily LZCep (HD 209481) is a double-lined spectroscopic bi- 
r-j ■ reach 60% (e.g. ISana et all l2008h . The effects of binarity on nal T s y stem composed of a primary component with a spec- 
k's ' stellar evolution depend not only on the properties of the com- ttaLtype O 8.5III and an O 9.5V secondary (|Conti & Alschuleij 
& ' ponents, but also on the mass ratio, the separation, the eccen- U2Z1I)- The projected rotational velocities were reported to be 
tricity. Consequently, their study is co mplex, requiring the ex- close to 155 km s and 105 km s by|Howarthet aL| <[1991|) . 
ploration of a large parameter space dde Mink et al.1 1201 lh . In The period of the system was found to be 3 070508 d by 
the case of young systems with large separations, the compo- iKameswara Rao ( 19_72|L and was confirmed by [Howarth et al.| 
nents of t he binary evolve independen tly, following single-star i l99l V and |Harries et al.| dl99J). The eccentricity estimated by 
evolution dPavlovski & Hensbergel2005h . In evolved binary svs- |Howarth et al.| dl99j is about 0.031. Although small, this is 



terns, the evolutionary scheme is completely different. If Roche significant enough to not be neglected. Indeed, according to 

lobe overflow happens, the surface chemistry and rotational rate |Lucy & Sweeney! 4l97li), the null hypothesis of a circular orbit 

of the components are deeply affected. In particular, the mass can be rejected with a significance level of about 0.05. However 

gainer is exp ected to show large nitrogen surface abundance en- lHarries et alj Q998J) determined an orbital eccentricity close to 

hancements (lLanger et al.l 120081) . In order to better understand zero in their orbital solution (e = 0.017 + 0.084). It thus appears 

the physics of binary systems and to test the predictions of the- that the circularization of the system is not clearly established, 

oretical models, it is thus important to quantitatively analyze the From the photometric point of view, LZ Cep is not an eclipsing 

system but ellipsoidal variations are clearly visible in the light 

Send offprint requests to: L. Mahy curve. The photometric analysis indicates that this binary is an 
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Fig. 1. NARVAL spectra between 4450 and 4570 A. The upper part of the figure shows the daily spectra. The lower two spectra are 
the results of the disentangling process. The two vertical dash lines correspond to the rest wavelengths of the Hei 4471 and Hen 
4542 lines 



evolved system dHill et al.lll976l) with a secondary component 
that fills up its Roche lobe. Hence, the system is most likely in 
a configuration of post case A mass transfer i.e., a system where 
the mass transfer happened w hile both componen ts are still in a 
hydrogen-burning core phase (Howart h et al.lll99lt) . 

LZ Cep is thus a target of choice for the study of the effects of 
binary interactions and their impacts on massive star evolution. 
The present paper is organized as follows. Sect. [2] presents the 
journal of observations and the data-reduction technique. Then, 
we devote the Sect. [3]to the determination of the orbital solution 
and the spectroscopic classification. Using the radial velocities 
obtained from the orbital solution, we apply a disentangling pro- 
gramme to separate the individual spectrum of both components. 
In Sect. |4] we apply the CMFGEN atmosphere code on the sep- 
arated spectra to determine the main stellar and wind parameters 
as well as the surface abundances of each star. These parameters 
provide us with information about the stars' evolutionary status. 
The last step in our investigation is reported in Sect. [5]and con- 
sists in the analysis of the Hipparcos light curve. It provides the 
system inclination and consequently the real masses of the stars. 
All the parameters and results determined in the present paper 
are discussed in Sect. [6] Finally, we present the conclusions in 
Sect. |7] 

2. Observations and data reduction 

LZCep was observed with the spectropolarimeter NARVAL 
mounted on the Telescope Bernard Lyot (TBL) located at the 
Pic du Midi observatory in the french Pyrenees. Table Q] gives 
the journal of observations. A total of ten spectra over a pe- 
riod of 14 d were obtained between July 20 th and August 3 rd 
2009 (Fig[TJ. This corresponds to more than 3 orbital periods of 



the system. The signal-to-noise ratios per 2.6 km s velocity 
bin are between 300 and 1300 depending on weather conditions 
and wavelength. A spectral resolution of 65000 is achieved. The 
data were automatic ally reduced with the Libre ESpRIT package 
(iDonati et alJ [T997) installed at TBL. An exhaustive descri ption 
of the reduction procedure is given in IDonati et al.l (1 19971) and 
we refer to this paper for information. 

3. Disentangling and orbital solution 

Previo us analyses by iHowarth et al. (199l]) and lHarries et al.l 
(1 19981) focused on the orbital elements of the system. In the 
present paper, we go one step further and study the properties of 
the individual components. The first step is to separate the spec- 
tra of each star. We applied an iterative disentangling method 
which consists in using alternatively the spectrum of one com- 
ponent (shifted by its radial velocity) and to subtract it from the 
observed spectra to determine the mean spectrum of the other 
component. The separated spectra, corrected by the brightness 
ratio (see below), are shown in Fig. Q] (the two bottom spec- 
tra) for the range 4450^4570 A and in Figures [3] and [4] (black 
lines) for a larger spectral band. Thi s disentangling method, in- 
spired by Gon zalez & Levatol (120061) . requires the knowledge of 
approximative radial velocities as input. We first measured these 
radial velocities by fitting two Gaussian profiles to the observed 
spectra. The cross-correlation technique was subsequently used 
by the disentangling program to refine the radial velocities from 
the disentangled spectra. This method ensures that radial veloc- 
ities are obtained at all orbital phases, including those where the 
lines are heavily blended. The cross-correlation masks were built 
on a common basis of the Si IV 4089^4116, He I 4388, 4471, 
4713, 4921, 5015, 5876, 6678, CIV 5801-5812, III 5592 and 
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Table 1. Journal of the observations. 



Date 


HJD 


exposure time 


SNR 


RV P 


RV & 




[d] 


[s] 




[km s '] 


[km s" 1 ] 


20 Jul 2009 


2455033.48549 


2700 


330-750 


26.7 


-100.1 


24 Jul 2009 


2455037.48790 


3300 


630-1320 


-96.6 


211.8 


25 Jul 2009 


2455038.49492 


2700 


/ion i t\r\f~\ 

480-1000 


26.3 


-95.8 


26 Jul 2009 


2455039.49494 


2700 


430-900 


45.7 


-153.9 


27 Jul 2009 


2455040.49228 


2700 


160-350 


-98.4 


206.6 


28 Jul 2009 


2455041.48290 


2700 


550-1180 


10.0 


-66.0 


29 Jul 2009 


2455042.49901 


2700 


640-1330 


52.1 


-179.1 


30 Jul 2009 


2455043.47908 


2700 


560-1150 


-95.8 


199.2 


31 Jul 2009 


2455044.49895 


2700 


570-1230 


-4.4 


-40.3 


03 Aug 2009 


2455047.46335 


2700 


500-1080 


-23.8 


5.3 



Notes. The HJD is measured at mid exposure. The signal-to-noise ratio (SNR) is given in the wavelength range 4000-6700 A and depends on the 
exact wavelength. Radial velocities are obtained by cross-correlation and have an error of 5 km s . 



Ha lines. For the primary component, we added the Hell 4200, 
4542, 5412 lines indicati ng, by the way, its earlier spectral type 
dConti & Alschulerlll97ll) . The list of the radial velocities ob- 
tained by cross-correlation is reported in TableQ] 




Fig. 2. Radial velocities as a function of orbital phase. The solid 
(open) symbols are for the primary (secondary). The circles are 
the data of the pres ent paper whilst the squares are the data of 
lHarries et al. I (fl998l) . 



We checked these radial velocities by computing the orbital 
solution of LZ Cep. For this purpose, we took into accou nt the 
well-k nown period of the system (P or b = 3.070507 d, Petrie 
119621) and we used the Liege Orbital Solution Package (LOSP, 
Sana&Gosset, A&A, submitted). This software is based on the 
generalization of the SB1 method of Wolf e et al.l dl9 67) to the 
SB2 case along th e lines described in iRauw et aD, (l2000h and 
Sana et al .]d2006a). We assigned a weight of 1.0 to all the mea- 
surements. Figure [2] represents the fitted radial velocity curves 
as a function of the phase. We added, for compa rison, the ra- 
dial ve locities measured from the WE spectra by lHarries et al.1 
(1998). In order to compar e both datasets, it w as necessary to 
shift the radial velocities of lHarries et al.l ([1998) by -9.4 km s _1 



to obtain the same systemic velocity than that measured for our 
dataset. We represent these data by squares in Fig. [2] To is the 
phase of the conjunction with the primary star in front. We tried 
to include a non-zero eccentricity in the orbital solution but the 
fit residuals are smaller for a zero eccent ricity, indicating that 
the system is circul arized. Previous studies (Howarth et al. 1991; 
lHarries et al.l lT998) have shown that this system could be in a 
post case A mass transfer phase (see also Sect. [5). It is therefore 
likely that the tidal interactions have circularize d that system. 
Other differences between our solution and that of lHowarth et al.l 
(1991) are notably a smaller value of smaller minimum 
masses and a higher mass ratio. The orbital parameters are listed 
in the upper part of Table [2] 

Finally, in the frame of the study of wind parameters (see 
Sect. @), we also decided to separate the IUE spectra in the 
range 1200-1800 A. We thus determined, from the orbital so- 
lution, the values of the radial velocities corresponding to the 
IUE spectra. We fixed these radial velocities in the disentangling 
program to obtain the resulting spectra shown in Fig. (black 
line). Unfortunately, the disentangled spectrum of the secondary 
is too noisy to be used and to constrain the wind parameters. 

In order to pinpoint the nature of the components, we deter- 
mined their spectral types. We measured the equivalent widths 
by fitting Gaussian functions to the line profiles. We did that 
for the most deblended observed spectra and from the disen- 
tangled spectra. The spectral-type determination is based on 
the quantitative classification criteria for O -type stars proposed 
bv IConti & Alschulerl (fl97lt) . IContil (fl973h and iMathvsl (fl988l 
119891) . We derived, from the ratios between the appropriate lines, 
the spectral classifications 09III and ON 9.7V for the primary 
and the secondary star, respectively. The "N" classification of 
the secondary comes from the strength of the nitrogen lines - e.g. 
Nn 4037^11-43, Nil 4601-07-1 3 and N m 4523 -30-34. These 
results confirm the conclusion of Hilditch ( 1974) who reported 
that the secondary is of different luminosity class than the pri- 
mary component. Otherwise, the spectral types ar e sligh t ly dif- 
ferent from those quoted in the literature. Indeed. iPetrid 01962) 
quoted LZ Cep as a binary system composed by two 09.5V 
stars whilst lConti & Alschulerl (1 197 lb determined a spectral type 
8.5111 for the primary and 09.5V for the secondary compo- 
nent. 

In a second step, we estimated the light ratio. For that 
purpose, we measured, for both stars, the equivalent widths 
of the Hel 4026, 4143, 4388, 4471, 5876, Hell 4200, 4542, 
Si IV 4089, OIII 5592 lines and we then computed the ratio 
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between these lines and the corresponding t heoretical av erage 
values reported in IConti & Alschulerl (1 19711) . IContil (1 19731) and 
IConti et ail (1 19771) . Due to the absence of O 9.7 stars in the lat- 
ter three papers, we used an 09.5 star as reference for the sec- 
ondary. We find that the primary star is about 1.85 times brighter 
than the secondary (/1//2 = 1-85 + 0.1). Alternatively we also 
measured the equivalent widths of the above-quoted lines in syn- 
thetic spectra computed with a similar effective temperature and 
logg to those of the primary and secondary components esti- 
mated (see Sect. HJ. With these synthetic spectra, we obtained 
a light ratio of about hlh = 1.91 which is consistent with our 
previous value (hlh = 1.85 ± 0.1). Consequently, we ad opt in 
the fo llowing a light ratio of about hlh = 1.85 + 0.1. Petrie 
( 1962) reported a value of hlh - 0.66 + 0.1 which corresponds 
to hlh = 1.52 + 0.3 (the method they used is not describe d 
though). Another estimate was made by lHarries et al.l (Q998). 
Although their method appears relevant, the result seems ques- 
tionable, probably due to the confusion between the primary and 
secondary stars. Indeed, the authors reported a value of about 
hlh - 0.65 ±0.1. How ever, when w e comp are the equiva- 
lent widths computed by lHarries et al.l (1 19981) with our mea- 
surements, no significant difference appears, indicating that the 
equivalent widths of the primary are larger than those measured 
on the secondary. W e thus assume that the real determination of 
lHarries et al.l (1 19981) is hlh — 0.65 ± 0.1 which corresponds to 
hlh = 1-54 ± 0.3. Cons equently, the discre pancy betwe en our 
measurement and that of lHarries et al.l d!998l) and lPetrid d 19621) 
probably comes from the sample of spectral lines that is larger 
in our analysis, rendering more accurate the determination of 
the light ratio Q. As a comparison, if we use a smaller bright- 
ness ratio (hlh ~ 1-54), we obtain absolute luminosities of 
log(Li/L ) = 5.08 and log(L 2 /L e ) = 4.74. The effect on the 
determination of the abundances is however small and does not 
exceed the error bars quoted in Sect. [4] 

Adopting our brightness ratio, we computed the luminosity 
of each component. The magnitud e and color of the system ar e 
V = 5.55 and (B - V) = 0.06 (iMafz-Apellaniz et al.l 12004). 
The visual extinction directl y follows from (B - V )q = -0.27 
(appropriate for late O stars, Mar tins & Plezll2006h . With these 
photometric parameters and assumi ng a distance m odulus for 
LZCep equal to DM = 10.06+J}^ (ISota et al.ll20TTl and refer- 
ences therein), consistent with the value of 10.12 given by Petrie 
(119621) . we found a visual absolute magnitude of the entire bi- 
nary system to be equal to My — -5.53^ ^. We derived, by us- 
ing the brightness ratio, M Vp = -5.06+J^ and M Vs = -4.39+°^ 
and with the bolometric correctio ns corresponding to the spec- 
tral types of both components dMartins & Plez 2006), we ob- 
tained absolute luminosities of log(Li/L Q ) = 5.1 1^24? and 

lo gfLy/LtT)) = 4. 69t[} [5 (more details on the method are given 
in Mah^etaDHoiO). As a comparison, if we took into account 
a smaller brightness ratio (hlh ~ 1-54), we would obtain abso- 
lute luminosities of log(Li/L Q ) = 5.08 and log(L2/L ) = 4.74. 
Accordingly, the brightness ratio computed in the present anal- 
ysis does not affect significantly the determination of the results 
since these values remain within the measured error bars. 

4. Spectroscopic analysis 

The disentangled optical spectra of both components have 
been analyzed separately with the atmosphere code CMFGEN 

1 If we rely on the same spectral lines as Ha rries et al] fl.998) for the 
determination, we obtain a brightness ratio only slightly different from 
their estimate. 



Table 2. Orbital solution and stellar parameters. 





Primary 


Secondary 


P[d] 


3.070507 (fixed) 


e 


0.0 (fixed) 


T [HJD - 2450000] 


5032.019+0.002 


q\m\IM%) 


2.53+0.05 


y [km s ] 


-11.81+0.91 


-11.40+1.20 


A [KJll 3 J 


88.72+1.02 


224.48+2.58 


ci sin 1 [Rq] 


5.38+0.06 


13.61+0.16 


ivi sin 1 lMq\ 


7.00+0.18 


2.77+0.05 


rms [km s" 1 ] 


2.9198 


T e ff [K] 


32000+1000 


28000+1000 


1 L 


5 ir 019 


4 69 +0 19 


l°g£ 


3.5+0.1 


3.1+0.1 


M spec [M Q ] 


15.9?*; 


4 r 2 - 4 


M ev [Me] 


25 V 62 


18.0! 24 


radius [R Q ] 


1 1 7 +33 
ii.;_ 27 


g 4+2:6 
?.t 22 


He/H 


0.1+0.02 


0.4+0.1 


C/H [x 1CT 4 ] 


1.0+0.5 


0.3+0.2 


N/H [x 1(T 4 ] 


0.85+0.2 


12.0+2.0 


o/h [x itr 4 ] 


3.0+0.5 


0.5+0.3 


V sin; [km s -1 ] 


130+10 


80+10 


Vmac [km s _I ] 


40+5 


44+5 


M [1(T 8 M yr 1 ] 


1.0+0.3 




Voo [km s~'] 


1800+100 





Notes. The given errors correspond to 1-cr. The solar abundances for 
the chemical elements quoted here are He/H = 0.1, C/H = 2.45 x 10 -4 , 
N/H = 0.60 x 10~ 4 , O/H = 4.57 x 10~ 4 , respectively. 



(lHillier&Millerlll998l) . The UV spectrum of the primary was 
also used. CMFGEN produces non-LTE models including winds 
and line-blanketing. The rate equation are solved iteratively to- 
gether with the radiative transfer equation to provide the oc- 
cupation numbers of energy levels and the radiation field. The 
computations are done in spherical geometry and in the co- 
moving frame. The temperature structure results from the ra- 
diative equilibrium condition. The velocity structure is an input 
of the calculation. In our case, we used TLUSTY hydrostatic 
structures in the photosphere connected to a /3-velocity law for 
the w ind part. The OSTAR2 002 grid of TLUSTY models was 
used dLanz & Hubeny||2003l) . The density structure directly fol- 
lows from the mass-conservation equation. Once the atmosphere 
model is obtained, a formal solution of the radiative transfer 
equation is performed to yield the synthetic spectrum. For that, 
accurate line profile including Stark broadening and a depth vari- 
able microturbulence are used. The latter varies from 10 km s _1 
in the photosphere to 0. 1 xv M in the outer atmosphere (v M being 
the terminal velocity). The spectrum is subsequently compared 
to the observational data to constrain the main parameters. In 
practice, we proceed as follows: 

- effective temperature: T e ff is determined from the ionization 
balance of He. The relative strength of He 1 and He n directly 
depends on the temperature. We used Hei 4471 and Hen 
4542 as the main indicators. Other lines were also checked: 
He 1 4026, He 1 4389, He 1 4713, He 1 5876, He n 4200, He n 
4686, Hen 5412. An accuracy of ±1000 K in the r eff deter- 
mination was achieved. 

- gravity: larger log g broadens the width of the Balmer lines 
wings. We thus used HB, Hy and H6 as the main gravity 
diagnostics. The final logg are accurate to 0.1 dex. 

- luminosity: the knowledge of the effective temperature, of 
the absolute V magnitude and of the extinction provides a 
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direct estimate of the luminosity as described in Sect. [3] The 
uncertainty on the luminosity determination depends mainly 
on the uncertainty on the distance (see Sect. [3}. 

- surface abundances: The surface abundances were deter- 
mined from the strength of He, C, N and O lines. He lines 
of all ionization states are stronger when the helium content 
is higher. This effect is different from a change in r e ff which 
changes the relative strengths of lines from successive ion- 
ization states. We have thus used the same lines as for the 
effective temperature determination to constrain the He con- 
tent. For nitrogen, we mainly relied on several N n lines be- 
tween 4-000 and 4700 A as well as the N m lines in the 4500- 
4520 A range. The main carbon indicators were the C m lines 
around 4070 A and C in 5696 A. For oxygen, we relied on 
the numerous O n lines between 4000 and 4600 A as well as 
on O in 5592 A. The typical uncertainty on the CNO surface 
abundance is 50%. To estimate these uncertainties, we ran 
models with element contents bracketing the preferred value 
and compared the corresponding spectra to the observed line 
profiles. When a clear discrepancy was seen on all lines of 
the same element, we adopted the corresponding abundance 
as the maximum/minimum value of the possible chemical 
composition. An illustration is shown in Fig. [6] As a test, 
we conducted the same analysis for the primary with four 
additional sets of models with the following sets of param- 
eters: r eff = 31000 K and logg = 3.5; T eS = 33000 K and 
log g = 3.5; T eff = 32000 K and logg = 3.4; T eff = 32000 K 
and logg = 3.6. We then estimated the best N/H abundance 
for each set of models (including also the best fit combi- 
nation, r e ff = 32000 K and logg = 3.5) by means of a 
chi-square method. We then computed the average and 1— 
sigma dispersion of this set of five N/H values and found 
N/H = 8.5 + 0.18 x 10~ 5 . This is perfectly consistent with 
the value we estimate from Fig. [6] For the other elements 
and the secondary, we thus proceeded as described above, 
running models with various abundances for the best fit pa- 
rameters (T e ff, logg). 

- mass-loss rate: we used the UV resonance lines of the pri- 
mary to estimate M. C iv 1538-1545 and N v 1240 were the 
main diagnostics. The uncertainty on M is 0.1 dex. For the 
secondary, the UV spectrum was not good enough to allow 
such a determination but in order to compute the atmosphere 
model associated to this star, we took a M — 3 x 10~ 8 M 
yr _1 . We also used homogeneous models (i.e. no clumping 
was included), which provided a good fit of the UV lines. 
As a consequence, our determination should be seen as an 
upper limit. We also included X-ray emission so as to have 
a canonical value of log = -7 (Sana et al. 2006b). This 
turned out to be crucial to reproduce the N v 1240 profile. 

- terminal velocity: the blueward extension of the P-Cygni 
profiles of the primary UV spectrum are directly related to 
Voo- We thus used them to estimate this parameter. Voo is de- 
rived to within 100 km s _1 . 

- rotational and macro turbulent velocity: we used th e Fourier 
transform formalism dSimon-Dfaz & Herrerol2007D to derive 
the projected rotational velocity (V sin/). The macroturbulent 
velocity was obtained by convolving the theoretical line pro- 
files with a Gaussian profile to best reproduce the shape of 
He 1 4713. 

The best fits of both components are shown in Fig. [3] and 
Fig. |U The fit of the UV spectrum of the primary is shown in 
Fig- E] The derived parameters are gathered in Table [2] (bottom 
part). The quality of the fit is on average very good. Only a few 





Fig. 5. Best CMFGEN fit (red) of the disentangled UV spectrum 
(black) of the primary component of LZ Cep. 
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Fig. 6. Illustration of the N/H determination for the primary com- 
ponent of LZ Cep. Colored lines are models and the black line 
is the disentangled spectrum. We estimate that the best fit abun- 
dance has N/H slightly above 8.0xl0~ 5 and that models with 
N/H = 6.0 x 10" 5 and N/H = 1.0 x 10" 4 correspond to the 
extrema of the allowed N/H values. 



features are not reproduced. This is the case of the Si in lines 
around 4550 A. Increasing r e jf leads to a better fit, but in that 
case the other diagnostics are less well reproduced. This problem 
is frequently observed in our models, and the origin is not known 
at present. Hence, we decided not to rely on these lines in the 
present study. 
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Fig. 3. Best CMFGEN fit (red) of the disentangled optical spectrum (black) of the primary component of LZ Cep. No data are 
available in the immediate vicinity of He 1 4921 . 



5. Light curve 

To better understand the binary system LZ Cep, and in particu- 
lar, its geometry, we also performed an analysis of the Hipparcos 
light curve. This investigation was made with the NIGHTFALL 
programme^. This software is based on a generalized Wilson- 
Devinney method assuming a standard Roche geometry. We per- 
formed a fit of the light curve by minimizing the free parameters. 
We fixed the effective temperatures to 32000 K and 28000 K, for 
the primary and the secondary star, respectively, as determined 



2 for more details, see the Nightfall User Manual 

by IWichmand d!998h available at the URL: 
www . hs . uni -hamburg . de/DE/Ins/Per/Wichmann/Nightf all . html 



by the CMFGEN analysis (see Sect. 0J. The mass ratio and the 
orbital period are those obtained from the orbital solution (see 
Table |2J. We stress that the zero phase of the light curve fit cor- 
responds to the zero phase of the orbital solution. 

The light curve (Fig. [7} does not exhibit any eclipses but 
ellipsoidal variations are clearly visible. The depth of the sec- 
ondary minimum is v ery similar to that of the primary minimum. 
iHowarth et al.l d!991l) considered two different models to explain 
the geometry of this binary system: a first one where both stars 
fill their Roche lobe, and a second one with a semi-detached con- 
figuration. The emphasis of their analysis was put on the second 
model with the secondary star filling its Roche lobe. 



6 



L. Mahy et al.: The two components of the evolved massive binary LZ Cep 



1 

0.8 






i i 1 i i 
H6 


i i 1 i i 


i | i i 


i i 














0.6 


: 1 


1 LU 


i K ii i 


1 1 1 1 U 1 


1 U J 


U 1 : 


7Hel 


Hel Nil 


CIII SilV SilVHel Hel 
■ i 1 i i 


Hel CIII Hell SilV Nil Sill 
■ i 1 i i 


en on 

■ i i i 


OH HeH 



4000 



4100 



4200 



4300 



4400 



= 1 1 1 1 1 


1 1 1 1 

Hel 


| i i i i | 


1 1 1 


- l 1 1 1 1 1 1 1 — 


i i i 


1 1 = 






X 








T 




r ii i 


1 




U 1 1 1 


1 1 


OIl/NII 

1 1 1 IU 1 


i 




1. Nil Hel Oil 
- .... 1 




Mgll 


NIII NII/III Hell 


Silll 

■ i i 


Oil SilV NIII CIII 
, 1 .... 1 , 


Hell 
. . . 


Hel J 



X 
3 



■V 

CD 
N 

6 

s- 

o 
•z. 



1 

0.9 
0.8 
0.7 



4400 



4450 



4500 



4550 




■ i i i i 



Ul ' 

Nil Hel 

_i i i i I i i i i I i i i i I i i i i I i i i i L 



4600 



4650 



4700 



1.1 


=-i — i — r 


I | I I I I | 


i i i | i i i i | 


i — i — r-= 


1 

0.9 




Ha A 










" T 




0.8 










0.7 


- 


1 1 1 1 1 1 1 


i Hei i 

■ 1 1 1 1 1 1 1 1 


- 

■ 1 1 1 1 1 1 



4800 4850 4900 4950 5000 5050 



6500 6550 6600 6650 6700 



-i 1 1 1 1 1 1 r- 



-i 1 1 1 1 1 1 1 r- 



~ i 1 r 



1 

0.9 
0.8 



w VyW *'» 



Hell 



Nil 
_l_ 



LU I 

Nil CIII 
i , l_ 



u 

CIV 

J I 



5400 



5500 



5600 
X[A] 



5700 



5800 



5900 



Fig. 4. Same as for Fig.[3]but for the secondary component of LZCep. No data are available in the immediate vicinity of He 1 4921. 



In our analysis, the NIGHTFALL programme also gives 
us two situations but, at each time, converges towards a semi- 
detached system. The first possibility is a situation in which the 
primary component fills up its Roche lobe (Sol 1 in Table [3]). 
This configuration favors an inclination of about 50.1°. The pri- 
mary star would fill up its Roche lobe at about 98% whilst the 
secondary component at about 87%. The second possible con- 
figuration is a system where the secondary star fills up its Roche 
lobe (Sol 2 in TableO. This configuration is closer to that found 
in the literature. Under this assumption, we found an inclination 
of about 48.1°. The volumes of the Roche lobes are filling up at 
about 95% and 99% f or the primary and the secondary compo- 
nents, respectively. As iHowarfh et al.l d!991l) already mentioned 
and no matter the chosen configuration, the inclination is well 



between 45° and 55° but, with the present investigation, we note 
that we improved, with more accuracy, the true value of the incli- 
nation. In both configurations, we included the reflection effects 
which are not negligible when one star fills out a large fraction 
of its Roche lobe. We also note that there is no photometric evi- 
dence of a non-zero eccentricity. 

The error bars given in Table [3] are established by exploring 
the parameter space. We fixed one parameter and let the others 
vary to reach the minimum of the x 1 ■ The error bars were then 
determined for a variation of the^ 2 corresponding to a 90% con- 
fidence level and the appropriate nu mber of freedom degrees (for 
more details, see lLinder et al.ll2009l) . 

Both models are really close to each other in terms of param- 
eters values and it is difficult to favor one configuration without 
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Fig. 7. Top: Photometry of LZ Cep from Hipparcos satellite. The 
red triangles are the observational data. A period of ~3.07 d is 
used. The solid line corresponds to the theoretical light curve fit 
with the parameters listed in Table[3] Bottom: A 3D view of both 
stars when the secondary star fills up its Roche lobe. 

Table 3. Parameters fitted from the Hipparcos light curve. 



Parameters 


Sol 1 


Sol 2 


if] 


so. 


48 l +i u 


q{M l IM 2 ] 


2.53 (fixed) 


2.53 (fixed) 


Filling factor primary [%] 


98.0!!^ 


94.8! 2 ° 


Filling factor secondary [%] 


87.0^-0 


98.6+?j 


r c ff,p [K] 


32000 (fixed) 


32000 (fixed) 


Ta,s [K] 


28000 (fixed) 


28000 (fixed) 


M p [M ] 


15 5 +L0 


16 9 +1 ° 


M s [M ] 


6-1!!:!! 


6 7 +1 
"•'-is 


Rpole,\) [Ro] 


10.53J 


io.5 + ;- 2 


Rpole,S L^o] 


6 r 1 2 


7 l +i:2 




13 1 +L2 


12 4 +12 


Requ,S [Rq] 


6 9+1-2 


9 3 +12 



Notes. The index 'p' ('s') refers to the primary (secondary). R po \ e is the 
polar radius, and R eqll the equatorial radius. 



any eclipse and with only the light curve fit. We can only safely 
conclude that LZ Cep is in a transitory state where both compo- 
nents are close to entirely fill their Roche lobe. LZ Cep is close 
to be a contact binary. 

6. Stellar evolution and binarity 

The position of the two components of the LZ Cep system in the 
HR diagram is displayed in Fig. [8] The two stars are located sig- 
nificantly off the zero-age main sequence. The isochrones reveal 
that the primary would be assigned an age of 4-6 Myrs by using 
single star evolutionary tracks. The secondary would have an age 
of 7-9 Myrs. The distinct ages could be due to a different rota- 



tion rate. Indeed, as recently demonstrated bv lBrott et al.l (1201 ll) . 
isochrones based on rapidly rotating evolutionary models tend to 
be located closer to the ZAMS compared to isochrones based on 
tracks with lower rotation rates. The effect is present for stars 
with luminosities above 3.0 10 5 L in the Galaxy. But it is sig- 
nificant only for values of V sin/ larger than -250 km s _1 . Since 
the projected rotation rates we derive are 130 and 80 km s~'for 
the primary and secondary components, we can exclude rotation 
as an explanation to the age difference between both stars. 

We have determined various masses for the components 
of the LZCep system. From the orbital solution and the light 
curve analysis, we have derived the true present masses: about 
16 M and 6.5 M for the primary and the secondary compo- 
nents, respectively. Those masses can be compared to the spec- 
troscopic masses obtained from the luminosity, temperature and 
gravity. We obtained 15.9*?'! M Q for the primary and 4.1*|'^ 
M for the secondary. Within the error bars, the agreement is 
very good for the primary, and marginal for the secondary. This 
could be due to the deformation of the secondary (see Fig. |7]i 
and the fact that we derive only an average gravity. The evolu- 
tionary masses derived using single star evolutionary tracks are 
respectively 25.3*4? M e and 18. 0;^ M Q for the primary and 
the secondary. The former are in strong conflict with the lat- 
ter. Under normal circumstances, a star of initial mass around 
25 M Q does not loose 10 M th rough stellar winds befo re it 
enters the Wolf-Rayet phase (e.g. iMevnet & MaederlfeOOOl) . So 
normal evolution cannot explain the mass discrepancy for the 
primary. In addition, if the discrepancy concerned only one com- 
ponent, one could argue that mass was lost through Roche-lobe 
overflow and given to the other component. But in our case both 
stars have evolutionary masses much larger than present masses. 
This is an indirect indication that the system was affected by bi- 
nary interaction and that single-star evolutionary tracks cannot 
be used blindly to derive the initial masses of the components 
of LZ Cep. We also note that the spectroscopic masses we ob- 
tained ar e lower than what is expected from single-star evolution 
(see e.g. Martins et al. 2005). This is especially true for the sec- 
ondary, for which a mass of 4.1 + j | M Q is much lower than the 
10-15 M of iMartins et al.l {2005). Having such a low mass but 
still the luminosity and effective temperature of an O star imply 
the occurrence of processes not present in single-star evolution. 

Another indication of interaction within the system comes 
from the analysis of the surface abundances. In Fig.|9]we show 
the nitrogen over carbon ratio as a function of the nitrogen over 
oxygen. Again, sin gle- star evolutionary track s including rota- 
tional mixing from Mevnet & Maeder (120031) are shown. We 
see that the primary, located at the bottom-left corner of Fig. 
[9] displays surface chemistry consistent with single star evo- 
lution. However, the secondary features extreme nitrogen en- 
richment and carbon/oxygen depletion. Such high N/C and N/O 
ratios are only reached in the Wolf-Rayet phase in single star 
evolution (the Wolf-Rayet phase is indicated by bold symbols 
in Fig. |9). Interestingly, this is consistent with the high he- 
lium content of the secondary: He/H = 0.4, corresponding 
to Y = 0.6, is only observed in advanced states of evolution 
in single-star evolutionary tracks with initial masses of about 
20 Mq. This raises the question of whether the secondary is not 
in fact more evolved than the primary and is currently in a state 
of central helium burning. U sing the mass-luminosity relation 
of ISchaerer & Maederl fl992) for core H-free stars, we find that 
a star with a current mass of about 6 M (resp. 5 M Q ) should 
have a luminosity log j- = 4.85 (resp. 4.69). This is remark- 
ably consistent with the derived luminosity of the secondary 
(log £ = 4.69). 
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All indications thus favor the scenario in which the sec- 
ondary was initially more massive, evolved faster, transferred 
mass to the companion and is now close to be a core He-burning 
object. This could also explain the rather low gravity and mass of 
the secondary (if it is truly a main-sequence O star, log g should 
be much higher - see also above). However, there are two caveats 
to this scenario. First, the spectroscopic appearance of the sec- 
ondary is not that of an evolved massive star. With a spectral 
type ON 9.7V, it looks like a normal main sequence star. Said 
differently, despite of all the indications, the absence of strong 
mass loss characterizing evo lved massive stars (and affe cting 
spectral types) is challenging. IGrafener & Hamannl ((2008) have 
shown that the mass loss rate of Wolf-Rayet stars depended 
on the Eddington factor T e . Estimating this factor from the 
parameters we determined, we obtain T e ~ 0.17. This is 
rather low, and much lower than the range of values for which 
IGrafener & Hamannl (12008b provide mass loss prescriptions for 
Wolf-Rayet stars (their computations are restricted to V e > 0.4). 
Hence, the secondary may very well have a rather nor mal wind. 
If we use the mass loss recipe of IVink et al.1 d200lb . we ob- 
tain logM ~ -7.0, which is much more typical of O stars. 
Consequently, it does not seem impossible that the secondary 
is an evolved object with a rather normal wind and thus the ap- 
pearance of an O star. 

The second caveat is the chemical composition of the pri- 
mary. If almost all the envelope of the secondary was dumped 
onto the primary, the chemical pattern should be different. The 
primary should show evidence for N enrichment and CO de- 
pletion since material from the secondary would be processed, 
the secondary being significantly evolved. According to our de- 
terminations, there are only weak signatures of CNO process- 
ing. This favors a weak accretion of material from the donor. A 
way out of this puzzle is to invoke a very non-conservative mass 
transfer even though this kind of process was not obs erved in the 
two p ost Roche lob e overflow key systems : </>per dBozic et al.l 
119951) and RY Scuti dGrundstrom et al.ll2007l) . All in all, the sce- 
nario according to which the secondary was initially more mas- 
sive and evolved faster to the core-He burning phase while fill- 
ing its Roche lobe and transferring very inefficiently material 
onto the primary seems a valuable explanation of the system. 

But other explanation s can be invoked. fde Mink et al .1 (120071) 
and lLanger et al.l ([2008) have argued that strong nitrogen en- 
hancement can be obtained on interacting binaries. Two effects 
are at work: the mass gainer receives N-rich material usually 
from the inner layers of the donor; and the angular momentum 
received from the primary triggers rotational mixing which fur- 
ther increases the surface nitrogen content of the secondary. In 
total, at solar metallicities, the enrichment can be of a factor of 
ten compared to the initial N abundance. According to Table [2] 
the nitrogen surface abundance of the secondary star in LZ Cep 
is 25 times the initial value. This is qualitatively consistent with 
the predictions of binary evolution which depend on the orbital 
parameters and masses of the components. In that case, con- 
trary to the scenario described above, the primary is the mass 
donor. Since rotational mixing is triggered by mass transfer, one 
also expects helium enrichment in the secondary, which is what 
we ob serve in the LZCep. However, according to lde Mink et al.l 
(2009), to reach a mass fraction Y = 0.6 requires very short 
periods and occurs only on stars more massive than ~ 40 M Q . 
Another concern is that the rotational velocity of the secondary 
is not extreme. With V sin/ = 80 km s _1 , and assuming that the 
rotation axis is parallel to the orbital axis, this give s a ro tation 
rate of ~ 100 km s _1 . But as shown bv lLanger et al.l (120081) . the 
rotational velocity can reach very high values after mass trans- 




4.6 4.5 4.4 



log (T,„) 

Fig. 8. HR diagram with the position of the two components of 
the LZCep system shown by the red squares. The primary is 
the most luminous star. The evolutionary tracks computed for an 
initial rotational velocity of 3 00 km s _1 and the isochrones are 
from lMevnet & Maederl d2003l) . 



fer, and subsequently come back to more moderate values. The 
question is whether fast mixing can be triggered during this short 
spin-up phase. It seems unlikely, but further dedicated simula- 
tions should be run to tackle this problem. 

With the inclination, projected rotational velocities and ra- 
dius, we can estimate the rotational period of each star and check 
if synchronization has been achieved in the system. Assuming 
an inclination angle of 49.1° (the average of the two values 
from Table |3j, we obtain an orbital period of 3.09+0.24 d us- 
ing R = 10.50 R Q (average of the polar radii for the primary) and 
3.75+0.31 d using R = 12.75 R Q (equatorial radius). Similarly, 
the rotation period of the secondary is 3.15+0.45 d (polar ra- 
dius) and 4.21+0.60 d (equatorial radius). Compared to the or- 
bital period (3.07 d), the true rotational periods of both compo- 
nents (likely in between the polar and equatorial cases) are thus 
roughly 10 to 20% larger, implying that the synchronous co- 
rotation is not completely established in the system. However, 
given the small differences and the similarity of the primary and 
secondary rotational periods, the system is certainly on the verge 
of achieving it. 

All in all, we have thus gathered evidence for mass transfer 
and tidal interaction in the massive binary system LZ Cep. The 
scenario according to which the secondary was initially the more 
massive star appears to explain more observational constraints, 
although it is not free of challenges. What is clear is that single- 
star evolutionary tracks are certainly not relevant to explain the 
position of the two components in t he HR diagram. This was 
highlighted by Wellstein et al.l d2001l) who computed evolution- 
ary tracks for interacting binaries and showed that as soon as 
mass transfer occurs, dramatic deviations from single-star evo- 
lutionary track happen. Consequently, the evolutionary masses 
we derived in Table [2] are very uncertain. 
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N/O 



Fig. 9. N/C versus N/O for the components of LZ Cep (red 
squares). The inset zoo ms on the primary value . The evolu- 
tionary tracks are from iMevnet & Maederl d2003h . The dotted 
line (resp. dashed, dot-dashed, short-dash long-dashed lines) 
are tracks for a 15 M (resp. 20, 25, 40 M ) star. The bold 
li nes indicate the Wo l f-Ray et phase of the evolution according 
to lMevnet& Maederl (120031) . i.e. X(H) < 0.4 and T eS > 10000 
K. 



7. Conclusion 

We have presented a detailed analysis of the binary system 
LZ Cep. We used high-resolution and high signal-to-noise spec- 
tra collected with the instrument NARVAL on the TBL. We 
performed an orbital solution from which we did not find ev- 
idence for a non-zero eccentricity. We have applied a spectral 
disentangling programme to our data set. This gave us the in- 
dividual spectra of each component of the system. We have de- 
termined the spectral type of the system: 9III+ON9.7 V. We 
have subsequently performed a quantitative spectroscopic analy- 
sis with the CMFGEN atmosphere code. Although the secondary 
star presents a main-sequence luminosity class, its surface abun- 
dances show a strong helium and nitrogen enhancement as well 
as carbon and oxygen depletion. The primary barely shows any 
sign of chemical processing. The analysis of the Hipparcos pho- 
tometry has revealed ellipsoidal variations that can be explained 
by two configurations for this system: LZ Cep is a semi-detached 
system in which either the primary or the secondary almost en- 
tirely fills its Roche lobe. The ellipsoidal variations of the light 
curve correspond to an inclination of about 48°. Together with 
the spectroscopic orbital solution, this allowed us to determine 
the masses of the system: about 16.0 M for the primary and 6.5 
M for the secondary. These masses are in good agreement with 
the spectroscopic masses. The evolutionary masses obtained us- 
ing single-star evolutionary tracks are larger. Our results indi- 
cate that LZ Cep is in a transitory state following mass transfer 
from the secondary - which was initially the most massive com- 
ponent - to the primary. The secondary is now (close to be) a 
core He-burning object. The mass transfer must have been very 
inefficient to account for the weak chemical processing at the 
surface of the primary. The caveats to this scenario have been 
discussed. The system is not yet completely synchronized, but 



the individual rotational periods are only 10 to 20% larger than 
the orbital period. 

Our investigation of the binary system LZ Cep represents a 
significant improvement compared to previous studies. We have 
obtained the individual spectra of each component. This leads 
to the determination of a number of stellar parameters which 
could be used to constrain the evolutionary status of the system. 
This evolved system is thus particularly interesting since it pro- 
vides observational constraints to theoretical bina ry evolutionary 
models. Similar works exist for B stars (e.g., see lPavlovski et al] 
2009) but studies of systems composed of two evolved early- 
type stars are scarce. We have emphasized the need for theoret- 
ical binary system evolutionary models since comparison of the 
stellar properties of the components of evolved binary systems 
to single-star theoretical predictions revealed discrepancies. 
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